Pion Scalar Form Factor and the Sigma Meson 
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From recent analysis of the mr scattering amplitude, it has been claimed that there exists a broad 
and light a meson. However, if this meson really exists, it must also appear in other observables such 
as the pion scalar form factor. With the use of unitarity and dispersion relations together with chiral 
perturbation theory, this form factor is analyzed in the complex energy plane. The result agrees well 
with the empirical information in the elastic region and reveals a resonance pole at \fs = 445 — i235 
MeV. This gives further strong evidence for the existence of the a meson. 
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It is still controversial whether a broad and light isosin- 
glet scalar meson really exists. After having been omitted 
from the Particle Data Group for more than twenty years, 
this scalar meson has reappeared in the last two editions 
of the Particle Data Group under the entry /o(400-1200) 
or a pL This reintroduction is based upon several new 
theoretical analyses |§-f|], which have given support to 
the existence of a broad and light scalar meson. In ad- 
dition, there has also recently been some experimental 
indications of a production from central pp collisions || . 

The theoretical evidence for the existence of the a me- 
son comes mainly from model-dependent analysis of 7T7t 
scattering However, if this meson in fact exist, 

it must also appear in other processes containing 7T7t in 
the final state. Therefore, in order to further investigate 
whether the a meson really exists, other processes should 
also be analyzed. The pion scalar form factor is such a 
process which can give further important information on 
the a meson. This is somewhat similar to the well-known 
case of the p(770) meson, where the information on this 
resonance can be obtained both from tttt scattering and 
from the pion vector form factor. 

The pion scalar form factor has previously been cal- 
culated using the inverse amplitude method (IAM) HQ. 
This method is based upon the combination of unitarity 
and dispersion relations together with chiral perturba- 
tion theory (ChPT) and has been used in order to 
extend the range of applicability of ChPT. In particular, 
the IAM has been applied in order to account for possi- 
ble resonances since this method, contrary to ChPT, can 
produce resonance poles in the complex energy plane. In 
this brief report, the scalar form factor is investigated 
further using the IAM and the result is analyzed in the 
complex energy plane in order to find possible resonance 
poles corresponding to the a meson. 

The pion scalar form factor F is given by the matrix 
element of the quark density 



lmF(s)=a(s)F*(s)t° (s), 



(2) 



where cr(s) is the phase-space factor and ig is the isos- 
inglet scalar 7T7t partial wave. This relation implies that 
the phase of F will coincide with the 7T7T phase shift <5§ in 
accordance with Watson's final-state theorem p0| . The 
first important inelastic effect starts at around 1 GeV 
and is due to the KK intermediate state. Since the main 
interest will be in energies well below this inelastic effect, 
in the following the form factor will only be calculated 
using the elastic approximation. 

The scalar form factor has been calculated to two loops 
in ChPT both by a dispersive analysis JOl and more re- 
cently by a full field theory calculation |12|]. The result 
can be written as 



F(s) = F (0 \s) + + F (2 \s), 



(3) 



where F^ is the leading order result, F^ the one-loop 
correction, and F^ the additional two-loop correction. 
Since the form factor will be normalized to F(0) = 1 in 
the following, one has for the leading order term F^ = 1. 
The one-loop correction is given in terms of the single 
one-loop low-energy constant l\ ^| , whereas the two-loop 
correction contains the additional one-loop low-energy 
constants V[, Z£, and l\ together with the two-loop low- 



energy constants r r S2 and 



S3 



|12|. The superscript r 



indicates that these low-energy constants depend on the 
renormalization scale [A, whereas the full form factor is 
scale- independent. Since ChPT is a perturbative expan- 
sion, the unitarity relation (||) will only be satisfied per- 
turbatively 



ImF^O) 
ImF (1) (s) 
Imf( 2 '(s) 



= 0, 



= a(s) [rcF^ (s)t° (0) (s) + Ret° (1) (s)l . (4) 



(n i (p 2 )\uu + dd\n J r (pi)) = S^F(s), 



(1) 



where s — (p 2 — Pi) 2 - This form factor is analytical in 
the complex s plane with a unitarity cut starting at the 
7T7T threshold. In the elastic region the unitarity relation 
is given by 



This perturbative unitarity will restrict the applicabil- 
ity of ChPT to the very low-energy region. However, 
with the use of the IAM, the range of applicability of 
ChPT can be substantially extended to also include res- 
onance regions. The starting point for this method is 
to write down a dispersion relation for the inverse of 



1 



the form factor T — l/F MJtfj. In this dispersion rela- 
tion the unitarity relation (||) gives Imr = —ImF/\F\ 2 = 
—crt/F. Expanding this quantity to two loops in ChPT 
gives Imr = -a[t^(l - RcF W) + Re^ 1 )] which can 
be used on the unitarity cut in the dispersion relation 
for r. The subtraction constants may also be evalu- 
ated by expanding the function T to two-loop order as 
r (2) = i_ F (i) + ^(i) 2 _ir(2). TnuS; neglecting possible 
zeros in the form factor, one has the following dispersion 
relation 



1 



F(s) 



1 + ai-s + a 2 s 



ds'a(s') x 



tjfV) [l-RcFW(s')] +Re^ (1) ( S Q 
s' 3 (s' — s — ie) 



(5) 



where three subtractions are used in order to make the 
dispersion integral convergent. This relation can be sim- 
plified by writing a dispersion relation for the function 
I^ 2 ). Using perturbative unitarity §Q) one finds the this 
dispersion relation will be exactly similar to the one given 
in Eq. (|J). Thus, the IAM to two loops in the chiral ex- 
pansion gives the form factor as 



F(s) 



1 - FW(s) + F^r(s) - F( 2 )(s) 



(6) 



This expression for the form factor is formally equiva- 
lent to the [0,2] Pade approximant applied on ChPT and 
will therefore coincide with the chiral expansion up to 
two loops. However, with the IAM the range of applica- 
bility of ChPT is substantially extended. This is based 
upon the fact that the expansion oit/F used in the IAM 
works well over a much larger region than the correspond- 
ing expansion of F*t used in ChPT. In fact, the former 
expansion works well throughout the elastic region, even 
when the form factor has a resonant character 0. 

However, the IAM may generate poles on the physical 
sheet which violate the analyticity requirement. These 
poles are caused by the high-energy part of the disper- 
sion integral in Eq. (^|). Since this part is not expected to 
be well approximated, it may cause the right-hand side 
of Eq. (Q) to vanish and thereby generate spurious poles 
in the form factor. These poles should in principle be 
removed without any significant influence in the region 
of applicability of the IAM. A rather general method to 
remove possible poles and thereby restore analyticity is 
to put the imaginary part of the IAM back into a disper- 
sion relation. With three subtractions the result can be 
written as 



3 r°° 



F(s) = l + ±(r 2 )s + cs 2 + 



lmF(s')ds' 
7T J 4M 2 s' 3 (s' - s - ie) 



, (7) 



where both the subtraction constants and Imi* 1 are cal- 
culated from the IAM (||). Here, it is assumed that ImF 
does not contain any poles on the unitarity cut and three 
subtractions are used in order to suppress the high-energy 



part of the dispersion integral. Without any poles Eq. 
(0) is just an identity, but with poles the output will in 
general be different from the input. However, in the re- 
gion where the IAM is applicable, the difference between 
the form (]?]) and the form (|^) should be small. In fact, 
this method to remove possible poles is equivalent to the 
subtraction of the poles on the physical sheet from the 
original IAM This will be discussed in more detail 
elsewhere Q. 

The IAM to two loops depends on a number of low- 
energy constants which have to be determined phe- 
nomenologically. Unfortunately, the scalar form factor 
is not directly accessible to experiment. However, in the 
elastic region the phase of F is given by the itit I = 
S phase shift Sq , which is known experimentally. Fitting 
these phase shifts up to 0.9 GeV and using the 

value of the pion scalar radius (r 2 ) = 0.60 fm 2 fl]|J|, 
some of the low-energy constants in the IAM to two loops 
have previously been determined without taking possible 
poles into account [||. Here, in order to remove spurious 
poles, this fit is repeated with the form factor given by 
Eq. ([?]). The result is shown in Fig. |l| together with the 
experimental tttt phase shifts, from where it is observed 
that the IAM agrees rather well with the main bulk of 
the data all the way up to 0.9 GeV. Thus, the IAM sat- 
isfies Watson's final-state theorem [[[o| quite well in the 
whole elastic region. This fit gives the following values 
for the low-energy constants 



l\ = 1.53 x 10~ 3 , 
r r S2 = 2.25 x 10" 3 , 
r r S3 = 7.60 x 10~ 5 



(8) 



at the renormalization scale /i = M p — 770 MeV. Since 
the experimental data are not very consistent with each 
other, there has not been assigned any error bars on these 
low-energy constants. However, the obtained values of l\ 
and r£ 3 agree rather well with the recent determination 
of these low-energy constants using two-loop ChPT jl2| . 
As for r r S2 , this low-energy constant has so far only been 
estimated on the basis of the resonance saturation hy- 
pothesis |l9| with a result (l2| that is somewhat smaller 
than the value obtained above. In the future, it might be 
possible to determine this two-loop low-energy constant 
from independent observables |2^] and thereby check the 
value obtained here. 

The scalar form factor can be defined in the whole com- 
plex s plane. Since it contains cuts starting at the irir 
threshold, this will involve different Riemann sheets. In 
the elastic approximation there are two Riemann sheets, 
which are defined according to the sign of the center of 
mass momenta q = \/s — 4M%/2. The first or physical 
sheet has positive values of Img, whereas the second or 
unphysical sheet has negative values of Iraq. The form 
factor given by either Eq. (||) or Eq. (^) can indeed 
be extended analytically to the whole complex s plane. 
This analytic continuation will involve infinitely many 
Riemann sheets since the cut in the IAM comes from 
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logarithmic functions. However, only two of these sheets 
correspond to the first and second Riemann sheet that 
the form factor should reproduce. 

In Fig. H the absolute square of the form factor (Q) is 
shown in the complex energy plane on the first Riemann 
sheet. On the real axis the result agrees very well with 
the result of a dispersive analysis Jll],[l8| , where the scalar 
form factor has been determined from the experimental 
tttt/KK phase shifts. Furthermore, the form factor (|t]) 
is analytic in the whole complex energy plane with the 
correct cut structure starting at the tttt threshold. This 
is different from the orig inal IAM (g) which generates a 
pole on the negative s axis. However, this pole is removed 
by using the form (^) without any significant influence on 
the result in the region shown in Fig. ||. 

From this figure it is also observed that around 0.4—0.5 
GeV the form of |F| 2 is somewhat reminiscent of a reso- 
nant structure. However, in order to investigate whether 
this form is really associated with a resonance, one has 
to consider the second Riemann sheet. On this sheet 
resonances are characterized by poles in the complex en- 
ergy plane, where the mass (Mr) and width (r_fj) of the 
resonance can be related to the position of the pole by 

y/S p ole = Mr - 1 — . (9) 

In Fig. H the absolute square of the form factor (Q) is 
shown in the complex energy plane on the second Rie- 
mann sheet. One finds that \F\ 2 indeed generates two 
complex conjugated poles corresponding to a broad and 
light resonance. In fact, the position of these poles is the 
same for the two expressions of the form factor given by 
Eq. (H) and Eq. (0), respectively. From the position of 
the pole, the mass and width of this a meson is given by 

M a = 445 MeV , T a = 470 MeV. (10) 

This compares rather well with the values M a — 470 
MeV and r CT = 500 MeV obtained in Ref. {§. However, 
rather different values for M a and have also been ob- 
tained with other theoretical models Therefore, it 
is important to reduce the model-dependence when the 
mass and width of the a meson are determined. The 
IAM is based solely on the use of unitarity and disper- 
sion relations together with ChPT. Therefore, within this 
approach, any model-dependence in the mass and width 
of the a meson is due to higher order terms in the chiral 
expansion together with the present uncertainties in the 
values of the low-energy constants. 

The IAM is in fact a systematic approach which can be 
applied to any given order in the chiral expansion. Origi- 
nally, this method was applied to the scalar form factor in 
the one-loop approximation H with a result that is for- 
mally equivalent to the [0,1] Pade approximant applied 
on ChPT. In this case the IAM depends on the single 
one-loop low-energy constant l\ which can be determined 
from the ratio F K / F r [ppC|] . With this low-energy con- 
stants fixed, the IAM to one loop agrees rather well with 



the empirical information up to about 0.5 GeV. The one- 
loop approximation also contains a pole on the negative 
s axis. However, this pole can be removed by using the 
method discussed previously without any significant in- 
fluence on the result in the elastic region. Extending this 
result to the whole complex energy plane, one finds that 
the IAM to one loop also generates a resonance pole, 
where the corresponding mass and width of this reso- 
nance are given by 

M a = 463 MeV , T a = 393 MeV. (11) 

Comparing these values with the values obtained from 
the IAM to two loops ( |Io| ) , it is observed that the masses 
are very similar, whereas the difference in the widths is 
somewhat larger. However, this difference is not signif- 
icant compared to the large uncertainty in the width of 
the cr meson given by the Particle Data Group ffl. In 
view of this the convergence of the IAM is satisfactory 
for both the mass and width of the a meson. Hence, the 
corrections to the values given in (^|) due to even higher 
orders in the chiral expansion are expected to be of little 
importance. There is also an uncertainty in the obtained 
values for M CT and T a coming from the uncertainties in 
the values of the low-energy constants. However, this ef- 
fect should also be rather small and could be estimated 
when the values of the low-energy constants in the IAM 
are determined more accurately. 

The IAM is not restricted to the scalar form factor but 
this method is in fact quite general and has been applied 
to other processes as well. In particular, the IAM has 
been applied to tttt scattering where it also generates a 
resonance pole corresponding to the a meson plfl . In 
fact, the mass and width of the a resonance obtained 
from tttt scattering, M a = 440 MeV and r CT = 490 MeV, 
agree very well with the values obtained in ([l0|). This 
strongly supports the consistency of the IAM and gives 
additional evidence for the existence of the a meson. 

To summarize, the pion scalar form factor has been cal- 
culated by the use of unitarity and dispersion relations 
together with the chiral expansion. In order to satisfy 
the analyticity requirement, possible poles on the physi- 
cal sheet are removed from this IAM. The result agrees 
well with both the experimental tttt phase shifts and a 
dispersive analysis in the whole elastic region. Making 
an analytic continuation of the scalar form factor to the 
complex energy plane, one finds a resonance pole corre- 
sponding to a broad and light scalar meson. The values 
for the mass and width of this a meson are obtained 
in a rather model-independent way, contrary to previous 
determinations where different theoretical models were 
applied. Indeed, any model-dependence in the obtained 
values for M a and should be rather small compared to 
the present uncertainty in these quantities. All this gives 
further strong evidence for the existence of the contro- 
versial broad and light a meson. 



3 



Electronic address : hannah@nordita.dk 
[1] Particle Data Group, C. Caso et ai, Eur. Phys. J. C 3, 1 (1998); Particle Data Group, R. M. Barnett et al, Phys. Rev. 
D 54, 1 (1996). 

[2] N. A. Tornqvist and M. Roos, Phys. Rev. Lett. 76, 1575 (1996). 

[3] S. Ishida et al, Prog. Theor. Phys. 98, 1005 (1997); S. Ishida et al, ibid. 95, 745 (1996). 

[4] R. Kaminski, L. Lesniak, and B. Loiseau, Phys. Lett. B 413, 130 (1997); R. Kamihski, L. Lesniak, and J.-P. Maillet, Phys. 

Rev. D 50, 3145 (1994). 
[5] GAMS Collaboration, D. Aide et al, Phys. Lett. B 397, 350 (1997). 
[6] T. Hannah, Phys. Rev. D 55, 5613 (1997). 
[7] T. N. Truong, Phys. Rev. Lett. 61, 2526 (1988). 
[8] S. Weinberg, Physica A 96, 327 (1979). 

[9] J. Gasser and H. Leutwyler, Ann. Phys. (N.Y.) 158, 142 (1984); Nucl. Phys. B250, 465 (1985). 
[10] K. M. Watson, Phys. Rev. 95, 228 (1954). 
[11] J. Gasser and U. G. Meissner, Nucl. Phys. B357, 90 (1991). 
[12] J. Bijnens, G. Colangelo, and P. Talavera, J. High Energy Phys. 5, 14 (1998). 
[13] T. Hannah and T. N. Truong (in preparation). 
[14] L. Rosselet et al, Phys. Rev. D 15, 574 (1977). 
[15] S. D. Protopopescu et al, Phys. Rev. D 7, 1279 (1973). 

[16] B. Hyams et al, Nucl. Phys. B64, 134 (1973); W. Ochs, Ph.D. thesis, Ludwig-Maximilians-Universitat, 1973. 

[17] P. Estabrooks and A. D. Martin, Nucl. Phys. B79, 301 (1974). 

[18] J. F. Donoghue, J. Gasser, and H. Leutwyler, Nucl. Phys. B343, 341 (1990). 

[19] G. Ecker et al, Nucl. Phys. B321, 311 (1989). 

[20] T. Hannah, Phys. Rev. D 54, 4648 (1996); J. Gasser and H. Leutwyler, Nucl. Phys. B250, 517 (1985). 
[21] A. Dobado and J. R. Pelaez, Phys. Rev. D 56, 3057 (1997). 



4 




Vs (GeV) 



FIG. 1. The phase 5q of the scalar form factor in the elastic region below yfs = 0.9 GeV. The experimental tvtt phase shifts 
are from Ref. (circles), Ref. ^| (squares), Ref. (rJ (diamonds), and Ref. |n| (triangles). 




FIG. 2. The scalar form factor \F\ 2 in the complex energy plane on the first Riemann sheet. The complex energy E is defined 
according to s = (ReE + ilmE) 2 . 
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FIG. 3. The scalar form factor \F\ 2 in the complex energy plane on the second Riemann sheet. 
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